Stochastic Process

HW9

如有错误,烦请指正.

1

image-20231115153658945

  1. 错误. 若随机过程 {N(t)} 是强度为 λ 的 Poisson 过程,则 EN(t)=λtt 有关不为常数,不符合平稳过程的条件.

  2. 正确.

  3. 正确. 若 Markov 链是有限状态,则它一定有正常返态.

    image-20231115154511131

    又因为其不存在平稳分布的充分必要条件是不存在正常返状态的不可约闭集,该 Markov 链一定存在平稳分布.

    image-20231115154705409

  4. 错误. 只是有概率停留在这个状态. (可能有其它的吸收态,等等.)

  5. 正确. To be added.

2

image-20231115154912215

见 HW3.

第一问:所有车辆的达到过程为强度为 9 的 Poisson 过程,则这个过程的第一次事件的等待时间 W1Exp(9),EW1=1/9. 注意到单位时间为 10 分钟,则第一辆车的到达的平均时间为 10/9 分钟.

第二问:记红、银、白车的达到过程分别为 N1(t),N2(t),N3(t),非白车的到达过程为 N4(t),则 N3(t) 是强度为 λ3=5 的 Poisson 过程,N4(t) 是强度为 λ4=λ1+λ2=4 的 Poisson 过程.

3

image-20231115164450751

有理谱密度:由谱密度是 ω 的非负实值偶函数知,S(ω) 形如

s0ω2n+a2n2ω2n2++a2ω2+a0ω2m+b2m2ω2m2++b2ω2+b0,

式中 s00. 又由于 R(0)>0,故 S(ω) 应在 [0,+) 上可积,从而 S(ω) 的分母不能有实根,分母多项式次数至少应比分子高 2 以及 s0>0.

4

image-20231115165846212

为避免混淆,记实平稳过程 {X(t)} 的自相关函数为

rX(τ)=25+41+6τ2,

而其协方差函数记为 RX(τ).

考虑自相关函数

rX(τ)=E[X(t)X(t+τ)]

关于 τ 的极限. 因为 τX(t)X(t+τ) 相互独立,所以

limτrX(τ)=limτE[X(t)X(t+τ)]=limτEX(t)EX(t+τ).

注意到 {X(t)} 为平稳过程,其期望为常数,记为 μ,则

limτrX(τ)=limτEX(t)EX(t+τ)=μ2.

又因为 rX(τ)25,τ,得 μ2=25,μ=±5.

再考虑 {X(t)} 的协方差函数

RX(τ)=E[(X(t)μ)(X(t+τ)μ)]=rX(τ)μ2=41+6τ2.

方差为

Var(X(t))=RX(0)=4.

5

image-20231115170826041

6

image-20231115171940380

image-20231115173039234

image-20231115173055774

记第 n(n1) 个顾客的服务时间为 Tn,则 TnExp(λ),且 ETn=1/λ=20,λ=0.05. 所以 N(t) 是强度为 λ=0.05 的 Poisson 过程.

(1)

记上午 8:00 为 t=0 时刻,则上午 12:00 时为 t=240.

因为

N(t)=N(t)N(0)P(λt),

所以

EN(t)=λt.

λ=0.05,t=240

EN(240)=0.05×240=12,

即上午 8:00 到 12:00 的平均服务顾客数为 12.

(2)

在上午 8:00 为 t=0 时刻的条件下,“顾客停留的时间”应当为“顾客的等待时间”. 记第 n 个顾客的等待时间为 Wn,则 WnΓ(n,λ),EWn=n/λ. 设长为 t 的时间内共服务了 N(t) 个顾客,则 N(t)P(λt),EN(t)=λt. 所求为

E[1N(t)i=1N(t)Wi]=E[E[1N(t)i=1N(t)Wi|N(t)]].

先算条件期望. 因为在给定 N(t)=n 的条件下,Wi,i=1,2,,n 的联合密度是与 (0,t] 上均匀分布中随机样本 Ui,i=1,2,,n 的次序统计量 U(i),i=1,2,,n 的联合密度是一样的,所以

E[1N(t)i=1N(t)Wi|N(t)=n]=E[1ni=1nWi|N(t)=n]=1nE[i=1nU(i)]=1nE[i=1nUi]=1nnEUi=t2.

所以(t=240

E[1N(t)i=1N(t)Wi]=E[E[1N(t)i=1N(t)Wi|N(t)]]=E[t2]=t2=120,

即这段时间内服务完的顾客停留的平均时间为 120 分钟.

image-20231115182223980

初始状态为 P0=(1/3,1/3,1/3).

P(n)=P0Pn.

(1)

第一个季度后三种牛奶的市场占有率向量为

P(1)=P0P=(1/3,1/3,1/3)(0.60.30.10.30.20.50.10.20.7)=(1/3,7/30,13/30).

第二个季度后(半年后)三种牛奶的市场占有率向量为

P(1)=P0P=(1/3,7/30,13/30)(0.60.30.10.30.20.50.10.20.7)=(47/150,7/30,34/75).

(2)

(3)

该 Markov 链是不可约的遍历链,则它存在唯一的极限分布,且此时的极限分布就是平稳分布,即

limnP(n)=π=(π1,π2,π3).

求解平稳分布:解方程组

{π1=0.6π1+0.3π2+0.1π3π2=0.3π1+0.2π2+0.2π3π3=0.1π1+0.5π2+0.7π3π1+π2+π3=1

π=(7/24,11/48,23/48).

image-20231115191506758

(1)

Xn1=i,则 Xn{1,2,,i}. 所以,

pij=P(Xn=j|Xn1=i)={0,j>i,1i,j=1,2,.i.
P=123N(10001/21/2001/31/31/301/N1/N1/N1/N).

(2)

N 个等价类,{1},{2},,{N}. 它们都是非周期的. 除 {1} 外的等价类都是非常返(瞬过)的,{1} 是正常返的.

(3)

状态 1 的平均返回时间

μ1=n=1nf11(n)=f11(1)+n=2nf11(n)=1+0=1.

因为状态 1 是非周期、正常返的,所以是遍历态,所以

limnpi1(n)=1μ1=1,i=1,2,,N.

因为状态 2,3,,N 是非周期的,所以

limnpij(n)=0,i=1,2,,N,j=2,3,,N.

所以极限分布存在,且为

limnP(n)=(1000100010001000).

image-20231115195027944

P(ω)=ω2+6,Q(ω)=ω4+8ω2+15=(ω2+3)(ω2+5)=(ω+5j)(ω+3j)(ω3j)(ω5j),Q(ω)=4ω3+16ω.

S(ω)=P(ω)Q(ω).

τ>0 时,由 Wiener-Khintchine 公式,

R(τ)=12π+S(ω)ejωτ dω.

借助留数定理,积分

+S(ω)ejωτ dω=2πjk=12Res[S(z)ejτz,ak],

其中 a1=3j,a2=5j 是复有理函数 S(z)=P(z)Q(z) 在上半平面内的全部极点,且都是一级极点.

对某个一级极点 ak

Res[S(z)ejτz,ak]=Res[P(z)ejτzQ(z),ak]=P(z)ejτzQ(z)|z=ak.

所以,

Res[S(z)ejτz,a1]=(z2+6)ejτz4z3+16z|z=3j=3e3τ43j=3e3τ4j,Res[S(z)ejτz,a2]=(z2+6)ejτz4z3+16z|z=5j=e5τ45j.

所以,

R(τ)=12π+S(ω)ejωτ dω=j[Res[S(z)ejτz,a1]+Res[S(z)ejτz,a1]]=34e3τ520e5τ,τ>0.

因为 R(τ) 是关于 τ 的偶函数,所以

R(τ)=34e3|τ|520e5|τ|.

image-20231115201810058

(1)

AU(1,1),EA=0,Var(A)=1/3,EA2=Var(A)+E2A=1/3.

X(t) 的均值函数

μX(t)=EX(t)=E[Asin(t+Φ)]=EAEsin(t+Φ)=0

为常数;

协方差函数

R(t,t+s)=E[(X(t)EX(t))(X(t+s)EX(t+s))]=E[X(t)X(t+s)]=E[A2sin(t+Φ)sin(t+s+Φ)]=EA2Esin(t+Φ)sin(t+s+Φ)=13E[12(cosscos(2t+s+2Φ))]=16(cossEcos(2t+s+2Φ))=16coss

只与时间差 s 有关,记为

R(τ)=16cosτ;

方差

Var(X(t))=R(0)=16,

二阶矩

EX2(t)=Var(X)+E2X(t)=16<

存在.

所以 {X(t),<t<+} 是平稳过程.

P1/21/2
Φπ/4π/4
2t+s+2Φ2t+s+π/22t+sπ/2
cos(2t+s+2Φ)sin(2t+s)sin(2t+s)
Ecos(2t+s+2Φ)=sin(2t+s)12+sin(2t+s)12=0.

(2)

考虑极限

limT1T02T(1τ2T)R(τ) dτ=16limT1T02T(1τ2T)cosτ dτ=16limT1T(02Tcosτ dτ12T02Tτcosτ dτ)=16limT1T(sin2Tsin2T+1cos2T2T)=16limT1cos2T2T2=0.

所以 {X(t),<t<+} 有均值遍历性.